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Abstract 



We solve a supersymmetric matrix model with a general potential. While matrix models 
usually describe surfaces, supersymmetry enforces a cancellation of bosonic and fermionic 
loops and only diagrams corresponding to so-called branched polymers survive. The 
eigenvalue distribution of the random matrices near the critical point is of a new kind. 



1 Introduction 



The theory of random matrices has been useful in the analysis of many different physi- 
cal systems. In particular in condensed-matter physics one encounters situations where 
physical quantities may not depend on the details of a (partly unknown) complicated 
Hamiltonian and where the average value of an observable can be calculated by replac- 
ing the average over eigenstates of the Hamiltonian with an average over an appropriate 
ensemble of random matrices. Aspects of vastly different physics such as electron local- 
isation phenomena in disordered conductors and semiconductors |l|, ^ H, disordered 
quantum wires and integer quantum Hall effects can be described by random ma- 
trix theory. Other applications can be found in nuclear physics where properties of highly 
excited nuclear levels can be described in random matrix language. Indeed, the very idea 
of random matrix theory started in nuclear physics with Wigner's seminal work and it 
was further developed by Dyson |^ and Mehta 0. For a recent review of the application 
of random matrix theory to nuclear physics see [|10|. More recently random matrices have 
been used in the study of the low energy chiral properties of QCD |jll| and to mention two 



extremes: the theory of random matrices has been used in the spectral theory of quan- 
tum mechanical systems with classically chaotic behaviour [jl2| and in two-dimensional 



quantum gravity. In the last case it has even played a derisive role in revealing a whole 
range of new non-perturbative phenomena (for recent reviews on this vast topic we refer 

to H). 

In many of the applications one is interested in computing ensemble averages which 
involve the density of states and correlation functions of densities of states of the Hamil- 
tonian H. The density is conveniently defined in terms of the trace of the resolvent of the 
Hamiltonian: 

G{z) = — ii—— (1.1) 
^ ' N z-H ^ ' 

where the normalization indicates that we have already replaced the real Hamiltonian 
with an X random matrix. The density of states for H is then given by 

p{E) = ^[G\E)-G-{E)l (1.2) 
where we have introduced the notation 

G^(E) = — tr i . (1.3) 

^ ' N E-H±ie ^ ^ 

It is now possible to calculate the average value over a number of microscopic energy 
levels by taking the average with respect to an appropriate ensemble of random matrices, 
typically a Gaussian ensemble of N x N Hermitian or orthogonal or complex matrices. 
For a general ensemble characterized by a potential V{H), one has for any observable 

ifiH)) = ^J dHfiH) exp[-NtrViH)l (1.4) 
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where the partition function Z is defined as 

Z = y"c//7 exp [-iVtrr(i/)]. (1.5) 

In [|l^ it was shown how to calculate Gi^z) in the large N limit for an arbitrary potential 
ViB.^ and the systematic 1/A^ expansion of Gi^z) was developed in [1^. In addition it is 
known that the large iV limit of the connected density-density correlator, or equivalently 



G{z,,z,) = U.- tr- -\ (1.6) 

\ Z\ — tl Z2 — n I conn 

is universal in the sense that it can be expressed as a function only of the endpoints of 
the distribution of eigenvalues {p{E)), independently of the form of the potential V. This 
was first shown in |]TB| and rediscovered in Essential aspects of the universality are 



valid to all orders in the 1/N expansion 

In many applications one needs both G^{E) and G^{E) and in these calculations it 
has been technically convenient to extend the random matrices to so-called supermatrices, 
which contain both Fermionic and Bosonic parts |jl8|, |l^ (for a review of the superma- 



trix theory as used in this context, see pOf). Even in the quantum gravity applications 



!1|, ^ as well as in purely combinatorial applications like the Meander problem 
supersymmetric matrix models have been most useful. 

In this paper we will show that a new kind of supersymmetric random matrix theory 
can be used to described a wide class of so-called branched polymers. 

Branched polymers are generalizations of the conventional linear polymers to systems 
where the polymer chain can branch at each node. The universal scaling behaviour of 
branched polymers is believed to be important for the understanding of macromolecules 
in solutions. Contrary to the linear polymers the scaling limit is not described by a 
conformal field theory In lattice statistical mechanics they are often modelled as 



"lattice animals", i.e. graphs of connected occupied sites on a lattice. Lattice animal 
exponents in d dimensions have been related to the Lee- Yang edge singularity m d — 2 
dimensions [^|. One usually introduces contact interactions heiween occupied neighbour 



sites which are not linked by a bond and solvent interactions between occupied sites and 
neighbouring empty sites. In addition one has a chemical potential associated with the 
growth of the branched polymers, i.e. the occupation of sites of the connected graph. 
One is interested in temperatures where the branched polymer becomes critical and the 
number of occupied sites grows to infinity. The detailed critical behaviour depends on the 
contact and solvent interactions as well as the entropy of the branched polymers. Part 
of the phase diagram can be approximated by taking the g — > 1 limit of an extension of 
the g-states Potts model |^ . 



The entropy factor is of primary importance for the phase transition. We use the word 
entropy in a generalized sense since it is intended not only to mean the number of branched 
polymers of a given size and shape, but also to take into account that a given branching 
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at a node will occur with a definite probability. In this paper we will concentrate on 
the analysis of abstract polymers, i.e. we ignore the embedding on a lattice and show 
how a supersymmetric matrix model of random matrices can be used to determine the 
universality classes of abstract branched polymers as a function of the weight attributed to 
branching at the individual nodes. This simplification implies that the partition function 
for our abstract branched polymer model can be defined as 

Z = ^u;(5P)e-^l^^l, (1.7) 

BP 

where \BP\ denotes the number of sites in the branched polymer and fi is the chemical 
potential. w{BP) is the weight of branching, i.e. it is a product of local weight factors 
associate with the branching at each site (or vertex): 

wiBP) = l[gis). (1.8) 

s 

We assume that g{s) only depends on the order of the site s, i.e. the number of neighbours 
connected to s by links and if s has order k we write for the weight. 

This model can be described by the following supersymmetric random matrix model: 

Z = J dWdW e-i''"^'"^'^'^\ (1.9) 

F(x)=x-^^a;^ (1.10) 
k>i ^ 

where W and W are "superfields" 

W, = {B,F), Wa = {B^,F), (1.11) 

a = 1, 2 and B and F are complex bosonic and fermionic (i.e. Grassmannian) N x N 
matrices, respectively. Supersymmetry means rotation between the B- and F-components 

6,B^ = Fe, 6,F = -eB, (1.12) 

SgB = eF, 6,F = -B^e, (1.13) 

where e and e are Grassmann valued matrices and it follows trivially that the model is 
invariant under the supersymmetric transformation ( |1.12|) - (|1.13|) . Note that it is a huge 
symmetry since e is matrix valued. 

The rest of this article is organised as follows: in section 2 and 3 we solve the random 
matrix model ( |1.9| ) using Schwinger-Dyson equations and supersymmetry Ward identities 
and show that the only diagrams which survive the supersymmetry are so-called cactus 
diagrams. In particular, we derive an equation which determines the possible universality 
classes. In section 4 we show how the cactus diagrams can be mapped onto branched 
polymer graphs and we give a purely combinatorial derivation of the important eq. ( p.l3|) . 



Section 5 discusses various generalizations of the model, while section 6 contains a short 
discussion of the results obtained. 
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2 Schwinger— Dyson equations and Ward identities 



The supersymmetric correlators of the model defined by the partition function (|1.9|) vanish 
due to the cancellations between bosonic and fermionic loops (see eq. (|2.5|)), but the 
correlators of the bosonic matrices are nontrivial. We shall study the generating function 
for these correlators: 

with 

Gn^^-triB^Bry (2.2) 
The Schwinger-Dyson equation for G{X) can be derived from the equality 

<'«'^»'^i:5|-(-^"*'"'«'"''«x3W).r« 



and reads 



^tTV'{WW)j^^'j=aX{{-tT^^] ). (2.3) 




To close this equation, we use the Ward identity which follows from the invariance 
of both the action and the measure in ( |1.9| ) under the supersymmetry transformations 



6, (f 1 1^bA) = 0, 

V u~WW \- B^B )l 

which, after a proper contraction of matrix indices, can be rewritten as 

WW 1 / 1 B^B \\ 
1 + — tr ■ 




WW X - BW \ N A - BW 



1 B^B 1 B^B \ , , 

77tr— — -). (2.4) 



WW X - B^B N A - B^B / ' 

This equation can be used to give a formal proof of the fact that the partition function 
is equal to unity and all supersymmetric correlators vanish. Taking the limit A ^ oo and 
retaining O (1/A) term in eq. (|2.4|) , one finds that 

(^tr [wwy^=0 (2.5) 

for all > 1. Since such correlators can be obtained by differentiating the partition 
function with respect to the coupling constants gk of the potential ( |1.10| ), the partition 
function is independent of the potential. 

On the other hand, the correlator entering the left hand side of eq. (^.3|) can be 



extracted from eq. ( p.4| ) in the large N limit, when factorisation holds. It can be seen 
expanding eq. (|2.4| ) in the powers of l/u. It is convenient to introduce the function 

Q(--^) = (w* %-U--A-flti; ) P,6) 
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and rewrite eqs. (|2.3|) and (|2.4| ) in the large limit as 

dui f dr] LjV'{rj)Q{rj,Lj) 



and 



aXG\X), (2.7) 



d7] 7]Q{r], f^^ f du_ ljG{uj) \ ^ r d^ ^Q{u,0 f dij_ ljG{uj) 



27ri V — 7] \ J 2Tii \ — uj ) J 27rz A — ^ J 27rz A — a; ' 

where the contours of integration encircle all singularities of Q and G, but not A, v and 
infinity. 

The integrals in eq. ( |2.8|) can be calculated by taking the residues at A, and oo. The 
residue at oo can be found since the asymptotic behaviour of Q{u,\) and G{\) follows 
from the definitions: 



and 



G(A) = i + o(i,). 
After some algebra one finds from 



This formula is to be substituted in eq. (|2.7|) . Then the integrals over 77 and uo can be 
done taking the residues at the poles of Q(?7, uj) as a function of r] and taking the residues 
aX oj = X and = 00. The result of this calculation reads 

(AG(A) - 1) V [\ - ^) = a\G\\). (2.10) 
It is now convenient to parametrise G{X) by 



and from ( 2.10 ) we obtain the following equation for x(A): 



xV\x) = (2.12) 
X — X 



If V{x) is a polynomial of n-th order, eq. ( p.l2|) is algebraic of (n + l)-th order. 
It follows from ( gTTD and (|2l|) that 
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Figure 1: Typical Feynman diagrams for (^jjtiB^Bj in the case where the potential 
V{x) = x - \g2X^ - \gzx^. 

so that, taking the limit A — > cxd in eq. (|2.12| ), we find the closed equation 

GiV\Gi)=a (2.13) 

for Gi. The surprising fact that we have obtained a closed equation for the propagator Gi 
is a consequence of the cancellations between bosonic and fermionic loops. In fig. |I]we have 
shown some of the diagrams which survive the cancellation. For obvious reasons we name 
them "cactus diagrams". Note that the diagrams have an orientation: the cactus loops 
can only proliferate on the exterior of already existing loops. This is in contradistinction 
to related cactus diagrams one encounters in the large N limit of purely bosonic or purely 
fermionic vector models p^ . 



In addition the Schwinger-Dyson equations for G„ defined by eq. (|2.2| ) allow us to 
express it in terms of Gk with k < n. We give the explicit formulas for the quartic 



potential V{x) = x — ^gx^- 



G. = (2.14) 



2aGn-l + «Z]fc+i=n-l GkGl + g I]r=2 ( — 1)"* Y.ki + ...+k,=n+l,ki<n Gki ■ ■ G 



r^2^G^ • ^^-^^^ 

To study the analytical structure of a;(A), it is useful to consider the inverse function 

x'^V'ix) 

A = I ^ , 2.16 

xV'[x) — a 

which is depicted schematically in fig. |^. It has a pole at x = Gi, as follows from eq. ( p.l3| ), 
a maximum at zero and a minimum at some xq > Gi. Consequently, the function x{X) 
has two branch points at A = and \ = z = A(xo). There may also be other branch 
points on the unphysical sheets of the Riemann surface of a; (A). The branch cut from to 
z is to be identified with the support of the eigenvalue distribution of the positive definite 
Hermitian matrix B'^B, the eigenvalue density being given by the discontinuity of G{X) 
across this cut. 
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Figure 2: Function \{x). The interval (0, z) represents an eigenvalue support of B'^B. 



3 Critical behaviour and scaling limit 



Inspection of eq. (|2.13|) shows that the critical behaviour comes about when the equality 

Gy{G,) = a, (3.1) 

holds simultaneously with 

[Gy{G,)^ = ... = [GX{Gy"'-''> = 0. (3.2) 
Near the critical point Gi behaves as 

G'i:^a-r'/"^(«c-«)'/"^, (3.3) 



where (5 = (—1)™-+^ [GcV'{Gc)] /m\. In general m = 2, but the multi-critical points 

dGi 
da 



with arbitrarily large m can be reached by tuning the potential. The susceptibility x ~ ^'^^ 



at the critical point scales as 

X ~ (oc - a)~^ 

with 7 = 1 — ^, which coincides with the (multi) critical index of the branched polymers 
30[| . In the next section we will make the mapping on a class of branched polymers 



explicit. 

Let us turn to the behaviour of the eigenvalue distribution at the critical point. Dif- 



ferentiating eq. ( |2.16| ) and expanding around ttc, one finds that Xq, defined by A'(xo) = 
(fig. 0), scales as xq — Gc + (c^c — ay^^"^~^\ Substituting this result into eq. ( |2.16| ), we 
obtain that z ~ Gcadc^c — <y)~^- Thus the eigenvalue distribution exhibits a behaviour 
which is rather unusual for matrix models — at the critical point it covers the whole 
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positive real semi-axis. This shows that an appropriate variable in the scaling limit is 

X(ac — a) 
P = 7; ' 

and that G{X) has the following form near the critical point: 



G(A) . i 




acGcP 



(3.4) 



The coefficient in front of the third term in the square brackets is chosen to match with 

From ( p.4| ) one can extract the scaling behaviour of the multi-point correlators of B^B: 

^ -f3-^{a,G,r~\a, - a)-"+i+^ /„ + d^iG,- (3.5) 

and f{p) is a generating function for /„. For the quartic potential the scaling behaviour 
( p.5| ) can be verified directly by using eq. ( |2.15| ). The generating function f{p) can be 
found by expansion of eq. ( |2.1(j| ) or eq. ( ^.12[ ) near the critical point. A simple calculation 
leads to 

1 / iV'"^ 

/(P) = - 1-- • (3.6) 
P\ Pj 

4 Mapping on branched polymers 

The graphs which survive the supersymmetric cancellation have an interpretation as 
branched polymers. Given a potential V{z) with coupling constants (see eq. ( |1.10| )) we 
have seen that the calculation of (^B^B^ to leading order in amounts to a summa- 
tion over all graphs of the kind shown in fig. |I|. The figure can be viewed as a branched 
polymer in the following way: (^B^b'^ consists of a closed loop of links with one of the 
vertices marked. The other vertices are created by expanding the interaction 



exp 



--A^tr (V{WW) - WW) 



a 



in powers of Qk and performing the Gaussian integrals. If the coupling constant is gk we 
attach a "A; — 1-fold blob" with weight factor gt/d if the coupling appears as 



^9k . 
a k 



(ww)' 



This procedure can now be iterated and we get a complete set of cactus diagrams 
which is in one-one correspondence with "chiral" branched polymers in the following way: 
the starting line and the end line of a cactus loop join at a vertex. Cut the loop open such 
that only the stating line is attached to the vertex. This produces a branched polymer 
graph which we call "chiral" since branching only occurs at one side of the open line, 
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Figure 3: The branched polymers corresponding to the cactus graphs shown in fig. |I]. 




+ 



Figure 4: The combinatorial equations for the branched polymer (or cactus) diagrams. 
The two equations lead to eq. ( p. 131) . 



corresponding to the fact that the cactus loops only can be attached to the exterior of 
already existing loops. This is illustrated in fig. ^. We can now write the following 
combinatorial identity for the chiral branched polymers (see fig. 



Gi = a + a'-^-V'^G^) + 1) + - V\G{) + l)' + 



a 



from which we get eq. (|2.13|) . 



(4.1) 



5 Ising model on branched polymer 



The purpose of the present section is to study a supersymmetric matrix model which 
describes the Ising model on the branched polymer. It can be constructed by introducing 
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a pair of the "superfields" and whose components transform under supersymmetry 
transformations by ( |1.12| ), ( |1.13| ) with the same parameters e and e for and W"^. We 
choose the interaction potential in the form 

A^tr (^^g+WiW^WiW^ + ^g.W2W'^W2W^ 

where g± = ge^^, and attach the spin variables to the vertices of the polymer obtained 
from the Feynman diagram for (^jjti BlB^'j as described in sec. The two types of 
vertices, the ones corresponding to g+ and (?_, are associated with spins of the opposite 
directions. 



Bl^B\^__ = ^^6u6,,D^^ (5.1) 



The Ising model is recovered if the propagators are defined by 

1 

Gauss _/Y 

with 

O: - ( ) . (5.2) 

The constants (5 and h are to be identified with the inverse temperature and an external 
magnetic field, respectively. 

In what follows we shall consider the model without the external field, so that the 
couplings gj^ and g^ are equal to each other. The matrix integral corresponding to this 
model has the form 

^- 2 , _ _ _ _ _ ^ 

Z = Yi dW^dW Q-^^'^{^^~'^^-^9WP+WWP+W-^gWP-WWP-W) g-j 

where 

= '-^ (5.4) 

and the inverse to (B.2|) matrix D"^ is 



^ " 2sinh2/? ( -e-^ ) ^^'^^ 



so that the quadratic part of the exponent in eq. ( |5.3|) involves 

1 



trWD-'W = -— tr 

2sinh2/5 



e^ (WiW^ + W2W^) - e-f^ (WiW^ + W2W^)] . (5.6) 



This model is invariant under the supersymmetry transformations, so the bosonic and 
fermionic loops are mutually cancelled, the partition function is equal to unity, and only 
the diagrams of the type depicted in fig. |l] survive. 

The Schwinger-Dyson equations for the two-point correlator of bosonic matrices can 
be obtained from the equality 



{[dw^dw^i-^Y: ^ 



a=l ij «J 



g - Af tr ( WD-^W-^ gWP+WWP+W- \ gWP- WWP- W) j^-y 

'^3 
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and can be rewritten in the form 



^tr [{BW-')^B^\)-g{^tr [{B^P+)^B^WP+W 
-gUtr Ub^P-)^B^WP-W 



(5.7) 



Similar to sec. we utilise the Ward identities to close eq. 

{5,{fPbIb^bI)) = q. 

Multiplying this equality by {P^Y^{P^Yi^ we get after some algebra the following equa- 
tion 



/— tr 

\n 



{B^P^)aBmP^W 



N 



tr 



tr 



{b^p^b 



{5.t 



Substituting it into eq. (pTTP, one obtains in the large N limit the closed equation for the 
2x2 matrix 

(5.9) 



PJ = (itri?t5n. 



From ( |5.7| ) and ( |5.8| ) we find 

D-^V -gSp (P+V) P+D - ^ Sp (P-V) P'V = 1. 



(5.10) 



For the symmetry reasons, the dressed coupling constants entering eq. ( p.lU| ) should 
be equal to each other: 

(5.11) 



g Sp (P+D j =gSp [P'V) = ^SpV = g. 



Multiplying eq. (|5.10|) by D from the left, we obtain 

V - gov = D. 



(5.12) 



Solving eq. ( |5.12| ) for V, we find 

1 



'.f^ -2 sinh2/3 g 



1 -2e/3^ + 2sinh2/?^2 



e^ - 2sinh2/3^ 



(5.13) 



Substitution of the solution (|5.13| ) into the definition of g (eq. (|5.11|) ) gives the equation 
from which it should be determined: 

2sinh2/?^3 - 2e ^ + (25sinh2/3 + 1)^ - ^e'^ = 0. (5.14) 

The third order equation ( |5.14| ), together with eq. ( p.l3[ ), completely determines the two- 



point correlators of the model. Which one to choose from the three solutions of eq. ( |5.14| ) 
is dictated by matching with perturbation theory. 
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Figure 5: The phase diagram of the Ising model on branched polymer. 



The critical behaviour comes about when two roots of the equation ( [5 .141 ) collide. It 
happens when the following conditions hold simultaneously: 

2smh.2/3g^-2e'^g^^ + {2g,smh2/3+l)gc- 9ce^ = 0, 

6sinh2/?^2 _4g/3^^^ (^2c/,sinh2/3+ 1) -^,e^ = 0. (5.15) 

The critical curve in the (3-g plane is depicted in fig. ^. The critical coupling constant 
behaves in the low temperature limit, P oo, as gdP) ~ In this limit the 

mixing between and W"^ is exponentially small and the above result reproduces the 
critical point of the solution (|2.14|) , ( p.l5| ) of the one-matrix model. The two-point 
correlators always have the square root singularity along the line of the phase transition; 
this recovers, in particular, the known result that the Ising model on the branched polymer 
never becomes critical |pl] . 



6 Discussion 



We have shown that the simplest non-trivial supersymmetric matrix model produces 
abstract branched polymers for the general potential ( p..l(]| ). Certain aspects of these 
abstract branched polymers have been found earlier by purely combinatorial methods, 
but the matrix model method allows a detailed study of the spectral properties. We 
predict that the non-trivial phases of the abstract branched polymers, corresponding 
the higher critical points with m > 2, also exist in "real" lattice animals, i.e. branched 
polymers embedded in regular d-dimensional lattices. Since these branched polymers 
look increasing similar to ordinary polymers for large m (see [^) there might exist a 
whole sequence of universality classes of lattice animals, interpolating between ordinary 
polymers with internal (Hausdorff) dimension 1 and "conventional" branched polymers 
with internal (Hausdorff) dimension 2. 
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In field theory the use of matrix models has usually been associated with Riemann 
surfaces. The large limit of Hermitian or complex matrix models has a diagrammatic 
expansion as graphs of spherical topology and an 1/iV expansion is an expansion in di- 
agrams with the topology of higher genus surfaces. Adding Supersymmetry to a purely 
bosonic or purely fermionic theory makes it more rigid and reduce the number of degrees 
of freedom. Here we observe this in an extreme way: only the tiny subclass of graphs 
corresponding to branched polymers survives for the model (|1.9| ), ( |1.10| ) when we look at 
observables like tr — B^B). 

It is an interesting question whether or not supersymmetric matrix models of the type 
discussed in this paper can be associated with super-Riemann surfaces. As is already 
mentioned in the previous paragraph, the potential ( p..lO| ) is separately invariant under 
both ( p..l2| ) and ( |1.13|) . This is the reason only cactus-like diagrams rather than planar 
graphs survive for the potential ( |1.1CI| ) in the large A^ limit. We can reduce the sym- 
metry by modifying the potential in order for a wider class of graphs to survive. Such 
supersymmetric matrix models deserve further investigations. 
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